Stability of negative and positive trions in quantum wires 
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Binding energies of negative and positive trions (X+) in quantum wires are studied for 

strong quantum confinement of carriers which results in a numerical exactly solvable model. The 
relative electron and hole localization has a strong effect on the stability of trions. For equal hole and 
electron confinement, is more stable but a small imbalance of the particle localization towards 
a stronger hole localization e.g. due to its larger effective mass, leads to the interchange of X~ and 
recombination lines in the photoluminescent spectrum as was recently observed experimentally. 
In case of larger X~ stability, a magnetic field oriented parallel to the wire axis leads to a stronger 
increase of the X^ binding energy resulting in a crossing of the X^ and X~ lines. 

PACS numbers; 68.65.La,71.35.Pq 



I. INTRODUCTION 

Trions are charged exciton complexes formed when an 
electron or a hole is boundi to a neutral exciton {X). 
The binding energies of the complexes are very small in 
bulk, but they are substantially enhanced in structures of 
reduced dimensionality, i.e., in quantum wellsSiMiSiSiL& 
and quantum wires^iifl*ii 

Due to the larger effective mass of the hole, in bulk^^ 
as well as in strictly two-dimensional confinemenlj^ the 
binding energy of positive trions [X'^) is larger than 
the negative trion {X~) binding energy. However, in 
quantum wells the observed^ X~ and X^ binding ener- 
gies are nearly equal, which is explained^'^ by a stronger 
hole localization within the quantum well enhancing the 
hole-hole interaction. The magnetic field perpendicular 
to the plane of confinement enhances more strongly the 
X~ stability leading to a crossing of X~ and X+ binding 
energiesiiii^ For trions localized on a defect of the quan- 
tum well potential X~ can become more stable than X^ 
even without the presence of an external magnetic field^ 
The combined quantum well and defect confinement cre- 
ates a three-dimensional potential similar to a quantum 
dot. In quantum dots the localization-related hole- hole 
interaction enhancement leads to the interchange of the 
order of the X~ and X^ recombination lines in the pho- 
toluminescence (PL) spectrum already for quantum dot 
diameters as large as 24 donor Bohr radii>i& For smaller 
dots the X'^ line becomes even more energetioiSiii than 
the X line. In coupled dotsi^ii^ this effect leads to the 
ground-state dissociation of X+jiLi^ for which the holes 
in the ground state occupy different dots. 

The present work is motivated by a recent experimen- 
tal studjiii of positive and negative exciton trions in 
V-groove GaAs/AlGaAs quantum wires. The negative 
trion was found to be distinctly more stable than X^ 
(binding energies oi X~ and X^ were determined as 4.2 
and 2.9 meV, respectively). Here, we indicate that the 
observedii order of X~ and X+ energy lines may be a 



consequence of modifications of the interactions due to 
a stronger hole confinement. In a previous theoretical 
study^ of trions in quantum wires X~^ was found to be 
more stable than X^, which was obtained in the case of 
equal hole and electron confinement. A crossing of X~ 
and A"+ PL lines as functions of the wire width has previ- 
ously been obtained in a quantum Monte-Carlo studjiiS 
of a quantum wire with a square well confinement poten- 
tial. In this paper we focus on the effect due to different 
electron and hole localization leading to modifications 
of the effective inter-particle interactions. We study the 
correlations between electrons and holes and consider the 
effect of a magnetic field oriented parallel to the quantum 
wire. The study of the stability of the trions is performed 
as function of the electron and hole localization instead 
of dimensions of the wire. It has been demonstrated^^ 
that in realistic quantum wires with strong confinement 
the binding energy of neutral excitons is governed by a 
size dependent parameter independent of the shape and 
composition of the wire. 

The stronger hole localization results from its weak 
penetration into the barrier material due to its larger 
effective mass than the electron band mass. For the 
spillover of the electron wave function out of the quan- 
tum wire, recently observed in self-assembled InAs/InP 
quantum wires^SS the ratio of the electron to hole local- 
ization can in principle be arbitrarily large. However, 
in the following we show that even a small enhancement 
of the hole localization changes the order of the X^ and 
X"*" PL recombination lines. 

For the purpose of the present study we apply the sin- 
gle band model for the hole and consider a harmonic 
oscillator confinement potential in the directions perpen- 
dicular to the wire, referred to as "lateral" in the fol- 
lowing. The present model does not account for the in- 
terface between the wire and barrier materials, so the 
effective mass discontinuity and dielectric constant mis- 
match are neglected. These effects usually strengthen 
the electron-hole interaction and weaken the penetra- 



2 



tion of the wave functions into the barrier. They are 
however of a secondary importance for GaAs/AlGaAs,^^ 
InAs/InP(SS and CdTe/ZnTe^'^ quantum wires. Note 
that, the present modelhng is inappUcable to the free- 
standing quantum wires, where the image charge effect 
is extremely strong^ 

We assume that the lateral confinement is strong, so 
that only the lowest subband for the electron and hole 
is occupied. This assumption allows for a reduction of 
the Schrodinger equation to an effective two-dimensional 
form. Usually the solution of the trion eigenequations 
is very challenging and requires extensive variational 
calculation o^i^i^i^i'^i^^i^^i-^'' or application of the quantum 
Monte Carlo methodspiS*i^ The present problem is unique 
in the sense that it allows for an exact inclusion of the 
interparticle correlations. 

The paper is organized as follows: the next Section 
contains the theory, the results are given in Section III, 
the conclusion and summary are presented in Section IV. 



II. THEORY 

We adopt the donor units, i.e., donor Bohr radius 
ad — iireoeh? /mee^ for the unit of length and twice the 
donor Rydberg 2Rd = /meO^ as the unit of the energy, 
where rUf, is the band electron effective mass and e is the 
dielectric constant. In these units, the Hamiltonian for a 
single electron in a quantum wire with harmonic oscilla- 
tor lateral confinement has the form 

(1) 

with the lateral Hamiltonian 
2 [dx^„ 



h; 



where le is the length of the harmonic oscillator confine- 
ment for the electron. The ground-state wave function 
of the Hamiltonian (2) is '^e = exp[-(a;^ -I- ?/^)/2Zg]/;e0r 
with the energy eigenvalue Ee — l/ll- In the adopted 
single-band approximation the hole ground-state wave 
function (vPft,) of the lateral confinement has the form 
of ^'e but with Ih - the harmonic oscillator length for the 
hole instead of 1^, and the energy is = l/cZ^, where 
a = mh/me is the hole to electron effective mass ratio, 
or, in other words, the hole mass in the donor units. The 
negative trion Hamiltonian can be written as 

1 1 1 

where Vgih {re2h) is the distance between the first (sec- 
ond) electron and the hole and ri2 is the electron-electron 
distance. We assume that the lateral confinement is suf- 
ficiently large that the trion wave function can be effec- 
tively separated into a product 

XX-{Zel,Ze2,Zh), (4) 



i7_ — H^i + Hp2 + Hh 



(3) 



where x- is the negative trion wave function of the mo- 
tion along the wire. The Hamiltonian (3) integrated over 
the lateral degrees of freedom with the wave function 
produces the effective trion Hamiltonian: 
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+ V''\le]Zel - Ze2) - {leh] Z^l ~ Z,,) 
-V'^leh; Ze2 - Zh), 



(5) 



with leh = vPi 
potential242^ 



-IfJ/'i and the effective interaction 



V^l-^z) = (7r/2)i/2erfc(|z|/\/2/)exp(zV2^^)/^, (6) 

which is finite at the origin {V°^{1;0) = l/l) and ap- 
proaches the 1/z asymptotic at large z. Hamiltonian 
(5) is written with respect to the sum of the ground- 
state energies of noninteracting two electrons and one 
hole. Therefore, the absolute value of the (negative) 
energy of a bound state is interpreted as the energy 
needed to separate all the particles away from one an- 
other. Introducing the center-of-mass coordinate Z = 
(zei + Ze2 + cFZh) /(2 -f cr) One obtains H'^ = -jM'iz^ + 
i?™', where M = 2 + ct is the negative trion mass and 
W"^^ is the relative motion Hamiltonian 
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+ V''{le;Z,,l~Z,,2)-V''{leh;Zhl) 

-V''\U;zh2), (7) 

with the reduced mass of an electron-hole pair fi = 
(t/(1 + cr), and the coordinates of the relative electron- 
hole positions Zhi = Zh — z^i and Zh2 = Zh — Ze2- In these 
coordinates the inter-electron distance along the length 
of the wire is Z12 = \zhi — Zh2\- The wave function x- is 
separable into a product of the center of mass and relative 
wave function X-izei, Ze2, Zh) = xcM{Z)x{zhi,Zh2)- 

The corresponding relative Hamiltonian for the posi- 
tive trion has the following form: 
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dzhidzh2 



+ V'''{lh;Zhl~ Zh2)-V''\leh;Zhl) 
-V'''{leh;Zh2), 



(8) 



with Zhi, Zh2 standing here for the relative position co- 
ordinates of the first and second hole with respect to the 
electron position. The reference energy for the Hamil- 
tonian |(SJ) is the energy of the dissociated complex, i.e, 
2Eh + Ee ■ 

In the following we consider also the exciton for which 
the effective Hamiltonian written with respect to the en- 
ergy of a dissociated electron and hole pair reads 



H 



X 
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(9) 
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The lowest eigenvalue of this Hamiltonian is equal to mi- 
nus the exciton binding energy {—E^). On the other 
hand the difference between -~E§ and the eigenvalues of 
trion Hamiltonians 17I8|I is equal to the trion binding en- 
ergies {Eg , Eg^) with respect to dissociation into an 
exciton and a free electron (for X^) or a hole (for 
Trion binding energies are equal to the red-shift of the 
trion recombination lines with respect to the exciton line 
in the PL spectrum. 
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FIG. 1: (a) [color online] Contour plot of the interaction po- 
tential V = V°f(L;zhi - Zh2) - V'iL-Zhi) - V''\L-Zh2) as 
function of the interparticle distances with lateral confinement 
length L = = = 1. Distances and energies are given in 
donor units. The dashed-dotted line corresponds to = 0. 
(b) The interaction potential plotted for L = 1 along the lines 
Zh2 = —Zhi, Zh2 = —Zhi — 1.5, and Zh2 = —Zhi — 3 marked 
in (a) with (thick) solid, dotted and dashed lines respectively, 
as function of the interelectron {X~) or interhole [X^) dis- 
tance zi2 = Zhi — Zh2- Thin solid line shows the — 3/|2:i2| 
asymptotic. 

We solve the relative Hamiltonian eigenequations us- 
ing the imaginary time technique^^ on a two-dimensional 
grid with a finite-difference approach. We use 201 points 
in both Zhi and Zh2 directions. The size of the computa- 
tional box in both directions is chosen " self-consistently" 
to be 12 times larger than the average distance be- 
tween the particles of the same charge defined as < 



{zhl - Zh2) >^l'^- 

In the present calculations we assumed harmonic oscil- 
lator lateral confinement which allows us to simplify the 
problem considerably because of the availability of an- 
alytical formula^'^^ for the effective one-dimensional in- 
teraction. The confinement lengths and Ih parametrize 
the strength of the particle localization. Since the single- 
particle energies cancel in the calculation of the trion 
binding energies, the applicability of the present results 
is wider. In fact the present results can be used for any 
form of the lateral confinement (which does not even 
have to be cylindrically symmetric) as long as it produces 
the same effective interaction potential. For instance the 
electron-electron interaction potential for = 2.95 and 
6 nm is very well (i.e. with a precision better than 2%) 
reproduced for a GaAs quantum wire (me=0.067) with a 
circular square well confinement of depth 320 meV and 
diameters 9.6 nm and 22.8 nm respectively. For elliptical 
harmonic oscillator confinement with different oscillator 
lengths in x and y directions {Ix and Zj,, respectively) 
we cannot give a closed analytical formula for the effec- 
tive interaction potential. Nevertheless, we have found 
via a numerical integration that the interaction potential 
between two electrons in an elliptical wire can be surpris- 
ingly well reproduced by formula (6) for a circular wire 
with an effective I = {lx + ^y)/'^- The numerically calcu- 
lated deviation between the two potentials is not larger 
than 2% for any interelectron distance. The essential as- 
sumption of the present model therefore does not rely on 
the form of the lateral confinement but on its strength, 
which has to be large enough to prevent the Coulomb 
interactions from deforming the lateral wave functions. 
The applied assumption of the frozen lateral degrees of 
freedom for the electron and the hole is applicable for the 
exciton binding energy when Zg < o,d and Ih < ad- This 
condition guarantees that the length of the lateral con- 
finement of the carriers is smaller than the bulk exciton 
radius, and that the sum of the lateral confinement ener- 
gies for the electron and the hole are at least two times 
larger than the exciton binding energy in bulk. For trions 
the applied approximation is better justified and the con- 
ditions are less stringent because the trions have a larger 
size and have smaller binding energies than the exciton. 



III. RESULTS 

We start the presentation of our results by discussing 
the properties of trions in quantum wires with equal lat- 
eral confinement for the electron and the hole (subsection 
III. A), and then in subsection III.B we show the effect of 
different confinements for the stability oiX^ and A"^ tri- 
ons. Subsection III.C describes the effect of the magnetic 
field oriented parallel to the axis of the wire. 
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FIG. 2: Wave functions for negative (a,d,g,j) and positive (b,c,e,f,h,i,k,l) trions for le = Ih = L = 0.2 in Zhi and Zh2 coordinates 
(horizontal and vertical axis, respectively) for different values of the mass ratio a. Plots (c,f,i,l) show the wave functions of the 
excited X"*" state antisymmetric with respect to the interchange of the holes. The dashed line in (f,i,l) shows the node of the 
wave functions. Plot (c) corresponds to an unbound state, for other plots the computational box is larger than the fragment 
displayed and the states are bound. 



A. Identical electron and hole lateral confinement 



For equal electron and hole lateral confinement {le = 
Ifi = L) the electron-electron, the hole-hole, and the 
electron-hole interactions have the same form. The total 
interaction potential, identical for both types of trions, 

is plotted in Fig. 1 (a) as fimction of z^i and Zh2 for 
L = 1. The regions of positive (negative) potentials are 



plotted with red (blue) colors. Zero of the interaction 
potential is marked with a dash-dotted line. The inter- 
action potential is minimal along the lines Zhi = and 
Zh2 = at which one of the two electrons and the hole are 
in the same position (for X^) or the position of one of 
the holes coincides with the electron position (for X~^). 
The potential is maximal along the diagonal Zhi = Zh2 at 
which the two particles of the same charge arc localized 
in the same point along the wire length. Fig. 1(b) shows 
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FIG. 3: Electron-hole (solid lines), electron-electron, hole- 
hole (dashed lines) pair correlation functions plots for X~ 
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FIG. 4: Binding energies of the negative (dashed lines) and 
positive trion (solid lines) states for L = 0.2 as function of the 
mass ratio a. Higher solid curve corresponds to the X'^ state 
antisymmetric with respect to the interchange of electrons and 
holes, ie. it is the first excited state. Thin vertical lines show 
the values of a = 1, 1.98, 6.72 and 15.2. Dotted curve, referred 
to the right axis, shows the exciton ground-state eigenvalue. 
Energies and lengths are in donor units. 



the cross-sections of the interaction potential along the 
three straight lines in Fig. 1(a) as function of the inter- 
electron {X~) or interhole (X+) distance zi2. On the 
antidiagonal the interaction has the form of a triangular 
potential well [cf. solid line in Fig. 1(b)]. Moving along 
the antidiagonal is equivalent to interchange the position 
of the two particles of the same charge with fixed posi- 
tion of the third particle of the opposite charge. Along 
the paths plotted with dotted and dashed hues in Fig. 
1(a), which are shifted below the antidiagonal, the po- 
tential has the form of a double potential well [see Fig. 
1(b)] with a barrier near the diagonal resulting from the 
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FIG. 5: [color online] Binding energies of the trions as func- 
tions of the length of the lateral confinement. Lines for 
a — 1, 1.98, 6.72 and 15.2 plotted with black, blue, red, and 
green colors, respectively. For a = 1 binding energies of X~ 
and X^ are equal. The dashed lines are the energies of the 
excited X"*" states antisymmetric with respect to the inter- 
change of the holes {X~ does not possesses a bound excited 
state for a > 1). Lines for X~ at cr = 1.98 and 6.72 have been 
omitted for clarity - they are situated between the a — 1 line 
and a — 15.2 line for X~(see Fig. 01 . Inset: energy eigen- 
values for neutral exciton and charged trions for a — 15.2. 
Energies and lengths are given in donor units. 



repulsion of the equally charged particles. For large zi2 
the potential approaches — 3/zi2 asymptotically which is 
shown by the this solid line in Fig. 1(b). 

Contour plots of the wave function of the negative and 
positive trions calculated for different effective mass ra- 
tios are plotted in Fig. 2 for L = 0.2. For Zg = Ih the 
negative and positive trion relative Hamiltonians (7) and 
(8) differ only by the factor standing in front of the mixed 
derivative (l/cr for and 1 for X+). The ground-state 
wave function for cr = 1 is the same for both trions [cf. 
Figs. 2(a,b)]. The first excited state of X~^, which is an- 
tisymmetric with respect to the interchange of the holes 
is unbound. Its wave function calculated for the size of 
the computational box 100 x 100 (in donor Bohr radius 
units) is plotted in Fig. 2(c). For the unbound state no 
computational box is large enough (i.e. the wave func- 
tion vanishes only at the ends of the computational box). 
The wave function is nonzero only near both axis. One 
of the holes stays at the position of the electron and the 
other strives to be as far as possible from the other two 
particles. 

Results for a — 1.98 plotted in Figs. 2(d-f) cor- 
respond to CdTe material parameters {mh = 0.19mo, 
TOg — 0.096mo) with the donor units ad = 5.4 nm and 
2Rd — 27.6 meV. The probability density maximum for 
(7 = 1 is split into two extrema at the antidiagonal of 
the plots [cf. Figs. 2(a-b)]. For a > 1 these two extrema 
merge into a single one for X~ [see Fig. 2(d)] and for X~^ 
they become more distinctly separated [see Fig. 2(e)] and 
the excited state for X+ becomes bound [cf. Fig. 2(f)]. 
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FIG. 6: [color online] Contour plot of the interaction potential for (a) negative trion V = V^^ {le; Zhi — Zh2) — V^leh', Zhi) ~ 
V°^{leh', Zh2) and (b) positive trion V = V^^ {lh \ Zhi — Zhi) — V^^ {leh \ Zhi) — V^leh', Zh2) as function of the interparticle distances 
for the lateral confinement lengths U = 1 and Ih ~ 0.5. Distances and energies are given in donor units. 
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FIG. 7: [color online] Shifts of the trion recombination PL 
lines with respect to the PL exciton line as function of the 
hole confinement length ilh) for GaAs. Different values of the 
electron confinement length are plotted with different colors. 



Figs. 2(g-i) and 2(j-l) for a = 6.7 and cr = 15.2 
correspond to GaAs (m/i — 0.45too, THe — 0.067mo, 
2i?d = 11.9 meV, — 9.8 nm) and InAs [m-t^ — 0.41mo, 
me = 0.027mo, 2i?d = 3.2 meV and ad = 29.7 nm) ma- 
terial parameters, respectively. Increasing a has an op- 
posite effect on the X~ and wave functions. For the 
negative (positive) trion the local minimum along the di- 
agonal Zhi = Zh2 is less (more) pronounced. The wave 
function evolution with a is related to the tunnelling of 
the particles of the same charge via the potential barrier 
presented in Fig. 1(b). The electrons in with light ef- 
fective masses tunnel easily through the diagonal barrier 
due to the interelectron repulsion. On the other hand the 
diagonal barrier is effectively much larger for the heavy- 
mass holes which prevents its penetratation at large a 
which leads to the appearance of the characteristic max- 



ima elongated along the diagonal in Figs. 2(e,h,k). 

The correlation between the particles in the complexes 
is more clearly visible in the pair correlation functions 
plotted in Fig. 3. The electron-hole correlation function 
is calculated as feh{z) = J dzhidzh2\x{zhi, Zh2)\'^S{z - 
Zhi) and the hole-hole (for X~^) , and electron-electron 

(for X-) as fsanie{z) = ^ dzhldzh2\x{zhl, Zh2)\^ 5[z - 

{zhi — Zh2)]- The Coulomb hole in the hole-hole correla- 
tion in X'^ is much larger than for electrons in X~ but at 
the expense of slightly weaker electron-hole localization. 

The binding energies of the exciton and the trions for 
L = 0.2 are plotted as functions of a in Fig. 4. All the 
binding energies are increasing functions of a. In bulk 
the first excited state of the positive trion is antisym- 
metric with respect to the hole interchange,'^^ possess 
the P symmetry and is bound for a > 4.2. The criti- 
cal value of the mass ratio is much smaller for quasi ID 
confinementi^ Here, for L = 0.2 the excited X+ state is 
bound for cr > 1.2 (see Fig. 3). For quasi ID confinement 
the lowest excited state has the S symmetry with respect 
to the axis of the wire but is of odd spatial parity, i.e., 
it is antisymmetric with respect to simultaneous change 
of sign of all the z coordinates (see Fig. 2). The ground- 
state of Ar+ becomes degenerate with respect to the sym- 
metry of the wave function, i.e., the hole interchange, for 
large a for which tunnelling through the diagonal poten- 
tial barrier (cf. Fig. 1) disappears. At large a also the 
probability density of the excited AT"*" level becomes iden- 
tical to the ground-state probability density (cf. Fig. 2). 
X~ does not possesses a bound excited state for cr > 1. 

The inset to Fig. 5 shows the ground-state energy for 
cr — 15.2 as function of the lateral confinement length. 
In the L — limit the average interparticle distances de- 
crease to zero and the energies diverge to minus infinity. 
This is a consequence of the Coulomb interaction singu- 
larity in one dimension. ™ The main part of Fig. 5 shows 
the shifts of the trion PL lines with respect to the exci- 
ton line (calculated as the difference of the eigenvalues 
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presented in the inset) for different values of cr. It turns 
out that the binding energies have a power law depen- 
dence on L, i.e. i^"^, for the and ground state 
presented in this figure q changes from 0.83 (cr = 1) to 
0.91 {X+ for a = 15.2). 



B. Effect of different electron and hole lateral 
confinement 

Let us now consider the interaction potential for 
stronger hole confinement. Figs. 6(a) and (b) show the 
interaction potentials for = 1 as in Fig. 1 but for 
smaller Ih = 0.5. For both the negative [cf. Fig. 6(a)] 
and the positive trion [cf. Fig. 6(b)] the potential min- 
ima at Zhi = and Zh2 = become deeper with respect 
to the le = Ih case presented in Fig. 1. For X~ the 
electron-electron interaction (the diagonal potential bar- 
rier) is not affected by the change of lu [compare Fig. 1 
and Fig. 6(a)]. On the other hand the hole-hole repulsive 
interaction for X^ is strongly increased. 




/g [nm] 

FIG. 8: [color online] Difference of the positive and negative 
trion binding energies (in meV) as function of the electron 
and hole confinement lengths for GaAs material parameters. 
Blue (red) regions correspond to more stable negative (posi- 
tive) trion. Above the dashed-dotted line the negative trion is 
unbound. The green line corresponds to Eb{X~) = 4.2 meV 
and the yellow line to Eb{X'^) — 2.9 meV. 

The effect of the hole localization on the trion binding 
energies is plotted in Fig. 7 for GaAs material parame- 
ters and fixed values of the electron lateral confinement. 
Consistently with the results of section III. A for l,, — Ih 
the positive trion is more stable than the negative trion. 
A decrease of In below the value of le results in the in- 
terchange of the X~ and X+ energy hnesiSS This is due 
to the enhanced hole-hole interaction shown in Fig. 6 
(b). The negative trion binding energy is a monotonous 
function of the hole confinement length, the larger 1^ the 
smaller is the electron- hole interaction stabilizing X^ . 
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FIG. 9: Interparticle distances for the negative (solid lines) 
and positive trions (dashed lines) and for the exciton (dotted 
line) for GaAs material parameters and le ~ 2.95 nm. (e-e), 
(e-h), and (h-h) stand for the electron-electron, the electron- 
hole and the hole-hole distance. 

The situation is more complex for X'^, since with in- 
creasing Ih also the destabilizing hole-hole interaction 
decreases. As a consequence the positive trion binding 
energy possesses a maximum as function of Ih- 

The difference of the positive and negative trion bind- 
ing energies is plotted in Fig. 8. Both the trions are 
equally stable for Ih = 0.92Ze — 0.38 nm. For Ih larger 
(smaller) than 0.921^ — 0.38 nm X~^ is more (less) sta- 
ble than X^ . The maximum of the AT^ binding energy 
presented in Fig. 7 follows a path that is nearly lin- 
ear for le > 2 nm and is approximately parametrized 
by Ih ~ 1.62le — 1.98 nm. For the points at the left of 
the dash-dotted line the electron-hole attractive interac- 
tion stabilizing the trion is so weak with respect to the 
electron-electron repulsive interaction destabilizing the 
complex that the negative trion stops to be bound (see 
also the hue for X~ at le = 2.95 nm in Fig. 7). The 
absence of the negative trion binding requires a substan- 
tially weaker hole localization than the localization of 
the electron which is rather impossible to obtain in the 
presently produced quantum wires and would require the 
valence band offset between the wire and the matrix to 
be much smaller than the conduction band offset. More- 
over, the present modelling based on the assumption that 
the lateral wave functions are not affected by the inter- 
action is likely to fail since the hole wave function is very 
likely to become more localized due to the attraction by 
strongly confined electrons. 

The fit of the calculated X^ and X^ binding energies 
to the experimental data is obtained at the crossing of 
the green and yellow lines, i.e., for le — 2.95 nm and Ih — 
1.3 nm. The obtained fit corresponds to realistic values 
which give a general idea on the particle localization in 
the wire (the measurementsii were performed on a V- 
groove GaAs/AlGaAs quantum wire with a thickness of 
the GaAs crescent of 3 nm at the center). Obviously, a 
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more realistic model is required to extract details of the 
confinement from the experimental data. 

The dependence of the size of the trion, i.e., the in- 
terparticle distance as function of the hole confinement 
length, is shown in Fig. 9 for 1^ = 2.95 nm. The electron- 
hole distance for the trions and the exciton have been cal- 
culated as ■y/<z^j>, and \/<z^h^, respectively. The 
hole-hole distance for X+ and th e electron-electron dis- 
tance for X~ are determined as ^ < {zhi — >■ The 
size of the exciton increases as 1^ increases which is due to 
the reduced value of the electron-hole interaction. Much 
stronger dependence on Ih is observed for , which be- 
comes unbound for Ih > 7 nm [cf. Fig. 7]. The de- 
pendence of the X'^ size on the hole confinement is non- 
monotonous. The positively charged complex has the 
smallest size near Ih — S nm when it is the most strongly 
bound [cf. Fig. 7]. For Ih — le ^ 2.95 nm the order 
of the interparticle distances in the two complexes is the 
same as in two-dimensional quantum wells (compare Fig. 
4 of Ref. fi2i] ) . In spite of the fact that the probability 
of finding both holes in X"*" in the same position is much 
smaller than for electrons in X~ (cf. Fig. 3) the longer 
tail of the electron-electron correlation function results 
in a larger electron-electron distance than hole-hole dis- 
tance. 




B[T] 



length {Ic = 36.28/^/B nm VT). Since 1(3) for high 
magnetic fields decreases to Ic, the magnetic field tends 
to equalize the electron and hole localization. The bind- 
ing energy of the trions can be obtained following paths 
of {lh{B),le{B)) on Fig. (8). The modification of the 
binding energies of the trions by the magnetic field is 
presented in Fig. 10 for different oscillator lengths. In a 
magnetic field of 40 T, 1{B ~ Q) ~ 4.9 nm is decreased 
to 1{B = 40T) = 4.4 nm and 1{B = 0) = 6.86 nm to 
1{B = 40T) = 5.19 nm. For 1^ = Ih = 6.86 nm (see the 
red curves in Fig. 10) the magnetic field decreases the 
length of confinement exactly as in the case presented in 
Fig. 5. In the more realistic case of stronger hole confine- 
ment, i.e., for le = 6.86 nm and Ih = 4.9 nm (see the black 
lines in Fig. 10) the magnetic field increases the binding 
energy of the positive trion more strongly. This is be- 
cause the magnetic field more strongly affects the larger 
le value than the smaller hole localization length Ih which 
increases the electron-hole interaction more strongly than 
the repulsive hole- hole potential. For X~ the effect of the 
increased electron-hole interaction is nearly cancelled by 
the increase of the electron-electron potential value. Near 
35 T a crossing of the black lines is observed, which cor- 
responds to passing from the "blue" to the "red" region 
on the phase diagram of Fig. 8. This crossing is qualita- 
tively opposite to the one obtained for two-dimensional 
quantum wells, ^"^ in which a small magnetic field (around 
1 T) increased the X^ stability over the X+ binding 
energy. In quantum wells the magnetic-field related in- 
crease of the single-particle energy is smaller for electrons 
and holes bound in the trion complex than for the low- 
est Landau level in the final state of the free electron 
and hole after the trion dissociation. The crossing of the 
binding energies observed in quantum wirea^ is there- 
fore at least partially due to the stronger dependence of 
the electron lowest Landau level - note an almost linear 
magnetic field dependence of the trion binding energies 
in Fig. 2 of Ref. [^^]. In the present calculations the free 
electron and the free hole are strongly localized in the 
plane perpendicular to the field and the single-particle 
magnetic field effects cancel due to the assumption of 
the frozen-lateral degrees of freedom, so that the cross- 
ing is entirely due to the modified effective inter-particle 
interactions. 



FIG. 10: [color online] Magnetic field dependence of the trion 
binding energies in GaAs for different values of the electron 
and hole oscillator lengths. 



C. Magnetic field parallel to the wire 

In the present approach it is straightforward to in- 
clude a magnetic field oriented parallel to the axis of the 
wire. It simply scales down the electron and hole oscil- 
lator lengths according to the formula 1{B) = {l/l'^{0) + 
l/lf)-^/^, where Ic = ^/2h/eB is the magnetic field 



In the case of stronger electron confinement {1^ = 4.9 
nm, Ih = 6.86 nm - the blue lines in Fig. 10) the situation 
is just opposite, the X~ is less strongly bound and the 
magnetic field acts more strongly on the negative trion 
decreasing the binding energy difference according to the 
mechanism described above. However, for the electron 
confinement equal or stronger than the hole confinement, 
the magnetic field does not lead to crossing of the trion 
energy lines and X~^ is more stable for any magnetic field. 
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IV. SUMMARY AND CONCLUSIONS 

We studied the properties of the negative and posi- 
tive trions in quantum wires with strong lateral confine- 
ment using the approximation of the lowest subband oc- 
cupancy which allows for a numerically exact solution of 
the multi-particle Schrodinger equation. We investigated 
the relative stability of the positive and negative trions 
with respect to the dissociation into an exciton and a free 
carrier for different electron and hole confinement. We 
found that the order of the negative and positive trion PL 
lines is interchanged when the lateral confinement of the 
hole is stronger than the one for the electron. In a GaAs 
quantum wire with If, = 5 nm we predict that when 
is 20% smaller the positive and negative trion recombi- 
nation lines interchange. The change in the order is due 
to modification of the effective interactions in the trion 
complexes. The present results provide an explanation 



for the recently experimentally observed larger stability 
of the negative trion in quantum wires We predict that 
for larger X~ stability due to stronger hole confinement 
the magnetic field oriented parallel to the axis will tend 
to interchange the order of the X'^ and X~ energy lines. 
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Binding energies of negative and positive trions (X+) in quantum wires are studied for 

strong quantum confinement of carriers which results in a numerical exactly solvable model. The 
relative electron and hole localization has a strong effect on the stability of trions. For equal hole and 
electron confinement, is more stable but a small imbalance of the particle localization towards 
a stronger hole localization e.g. due to its larger effective mass, leads to the interchange of X~ and 
recombination lines in the photoluminescent spectrum as was recently observed experimentally. 
In case of larger X~ stability, a magnetic field oriented parallel to the wire axis leads to a stronger 
increase of the X^ binding energy resulting in a crossing of the X^ and X~ lines. 
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I. INTRODUCTION 

Trions are charged exciton complexes formed when an 
electron or a hole is boundi to a neutral exciton {X). 
The binding energies of the complexes are very small in 
bulk, but they are substantially enhanced in structures of 
reduced dimensionality, i.e., in quantum wellsSiMiSiSiL& 
and quantum wires^iifl*ii 

Due to the larger effective mass of the hole, in bulk^^ 
as well as in strictly two-dimensional confinemenlj^ the 
binding energy of positive trions [X'^) is larger than 
the negative trion {X~) binding energy. However, in 
quantum wells the observed^ X~ and X^ binding ener- 
gies are nearly equal, which is explained^'^ by a stronger 
hole localization within the quantum well enhancing the 
hole-hole interaction. The magnetic field perpendicular 
to the plane of confinement enhances more strongly the 
X~ stability leading to a crossing of X~ and X+ binding 
energiesiiii^ For trions localized on a defect of the quan- 
tum well potential X~ can become more stable than X^ 
even without the presence of an external magnetic field^ 
The combined quantum well and defect confinement cre- 
ates a three-dimensional potential similar to a quantum 
dot. In quantum dots the localization-related hole- hole 
interaction enhancement leads to the interchange of the 
order of the X~ and X^ recombination lines in the pho- 
toluminescence (PL) spectrum already for quantum dot 
diameters as large as 24 donor Bohr radii>i& For smaller 
dots the X'^ line becomes even more energetioiSiii than 
the X line. In coupled dotsi^ii^ this effect leads to the 
ground-state dissociation of X+jiLi^ for which the holes 
in the ground state occupy different dots. 

The present work is motivated by a recent experimen- 
tal studjiii of positive and negative exciton trions in 
V-groove GaAs/AlGaAs quantum wires. The negative 
trion was found to be distinctly more stable than X^ 
(binding energies oi X~ and X^ were determined as 4.2 
and 2.9 meV, respectively). Here, we indicate that the 
observedii order of X~ and X+ energy lines may be a 



consequence of modifications of the interactions due to 
a stronger hole confinement. In a previous theoretical 
study^ of trions in quantum wires X~^ was found to be 
more stable than X^, which was obtained in the case of 
equal hole and electron confinement. A crossing of X~ 
and A"+ PL lines as functions of the wire width has previ- 
ously been obtained in a quantum Monte-Carlo studjiiS 
of a quantum wire with a square well confinement poten- 
tial. In this paper we focus on the effect due to different 
electron and hole localization leading to modifications 
of the effective inter-particle interactions. We study the 
correlations between electrons and holes and consider the 
effect of a magnetic field oriented parallel to the quantum 
wire. The study of the stability of the trions is performed 
as function of the electron and hole localization instead 
of dimensions of the wire. It has been demonstrated^^ 
that in realistic quantum wires with strong confinement 
the binding energy of neutral excitons is governed by a 
size dependent parameter independent of the shape and 
composition of the wire. 

The stronger hole localization results from its weak 
penetration into the barrier material due to its larger 
effective mass than the electron band mass. For the 
spillover of the electron wave function out of the quan- 
tum wire, recently observed in self-assembled InAs/InP 
quantum wires^SS the ratio of the electron to hole local- 
ization can in principle be arbitrarily large. However, 
in the following we show that even a small enhancement 
of the hole localization changes the order of the X^ and 
X"*" PL recombination lines. 

For the purpose of the present study we apply the sin- 
gle band model for the hole and consider a harmonic 
oscillator confinement potential in the directions perpen- 
dicular to the wire, referred to as "lateral" in the fol- 
lowing. The present model does not account for the in- 
terface between the wire and barrier materials, so the 
effective mass discontinuity and dielectric constant mis- 
match are neglected. These effects usually strengthen 
the electron-hole interaction and weaken the penetra- 
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tion of the wave functions into the barrier. They are 
however of a secondary importance for GaAs/AlGaAs,^^ 
InAs/InP(SS and CdTe/ZnTe^'^ quantum wires. Note 
that, the present modelhng is inappUcable to the free- 
standing quantum wires, where the image charge effect 
is extremely strong^ 

We assume that the lateral confinement is strong, so 
that only the lowest subband for the electron and hole 
is occupied. This assumption allows for a reduction of 
the Schrodinger equation to an effective two-dimensional 
form. Usually the solution of the trion eigenequations 
is very challenging and requires extensive variational 
calculation o^i^i^i^i'^i^^i^^i-^'' or application of the quantum 
Monte Carlo methodspiS*i^ The present problem is unique 
in the sense that it allows for an exact inclusion of the 
interparticle correlations. 

The paper is organized as follows: the next Section 
contains the theory, the results are given in Section III, 
the conclusion and summary are presented in Section IV. 



II. THEORY 

We adopt the donor units, i.e., donor Bohr radius 
ad — iireoeh? /mee^ for the unit of length and twice the 
donor Rydberg 2Rd = /meO^ as the unit of the energy, 
where rUf, is the band electron effective mass and e is the 
dielectric constant. In these units, the Hamiltonian for a 
single electron in a quantum wire with harmonic oscilla- 
tor lateral confinement has the form 

(1) 

with the lateral Hamiltonian 
2 [dx^„ 



h; 



where le is the length of the harmonic oscillator confine- 
ment for the electron. The ground-state wave function 
of the Hamiltonian (2) is '^e = exp[-(a;^ -I- ?/^)/2Zg]/;e0r 
with the energy eigenvalue Ee — l/ll- In the adopted 
single-band approximation the hole ground-state wave 
function (vPft,) of the lateral confinement has the form 
of ^'e but with Ih - the harmonic oscillator length for the 
hole instead of 1^, and the energy is = l/cZ^, where 
a = mh/me is the hole to electron effective mass ratio, 
or, in other words, the hole mass in the donor units. The 
negative trion Hamiltonian can be written as 

1 1 1 

where Vgih {re2h) is the distance between the first (sec- 
ond) electron and the hole and ri2 is the electron-electron 
distance. We assume that the lateral confinement is suf- 
ficiently large that the trion wave function can be effec- 
tively separated into a product 

XX-{Zel,Ze2,Zh), (4) 



i7_ — H^i + Hp2 + Hh 



(3) 



where x- is the negative trion wave function of the mo- 
tion along the wire. The Hamiltonian (3) integrated over 
the lateral degrees of freedom with the wave function 
produces the effective trion Hamiltonian: 



= - 
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+ V''\le]Zel - Ze2) - {leh] Z^l ~ Z,,) 
-V'^leh; Ze2 - Zh), 



(5) 



with leh = vPi 
potential242^ 



-IfJ/'i and the effective interaction 



V^l-^z) = (7r/2)i/2erfc(|z|/\/2/)exp(zV2^^)/^, (6) 

which is finite at the origin {V°^{1;0) = l/l) and ap- 
proaches the 1/z asymptotic at large z. Hamiltonian 
(5) is written with respect to the sum of the ground- 
state energies of noninteracting two electrons and one 
hole. Therefore, the absolute value of the (negative) 
energy of a bound state is interpreted as the energy 
needed to separate all the particles away from one an- 
other. Introducing the center-of-mass coordinate Z = 
(zei + Ze2 + cFZh) /(2 -f cr) One obtains H'^ = -jM'iz^ + 
i?™', where M = 2 + ct is the negative trion mass and 
W"^^ is the relative motion Hamiltonian 
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+ V''{le;Z,,l~Z,,2)-V''{leh;Zhl) 

-V''\U;zh2), (7) 

with the reduced mass of an electron-hole pair fi = 
(t/(1 + cr), and the coordinates of the relative electron- 
hole positions Zhi = Zh — z^i and Zh2 = Zh — Ze2- In these 
coordinates the inter-electron distance along the length 
of the wire is Z12 = \zhi — Zh2\- The wave function x- is 
separable into a product of the center of mass and relative 
wave function X-izei, Ze2, Zh) = xcM{Z)x{zhi,Zh2)- 

The corresponding relative Hamiltonian for the posi- 
tive trion has the following form: 



™i ^ — L 
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+ V'''{lh;Zhl~ Zh2)-V''\leh;Zhl) 
-V'''{leh;Zh2), 



(8) 



with Zhi, Zh2 standing here for the relative position co- 
ordinates of the first and second hole with respect to the 
electron position. The reference energy for the Hamil- 
tonian |(SJ) is the energy of the dissociated complex, i.e, 
2Eh + Ee ■ 

In the following we consider also the exciton for which 
the effective Hamiltonian written with respect to the en- 
ergy of a dissociated electron and hole pair reads 



H 



X 
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(9) 
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The lowest eigenvalue of this Hamiltonian is equal to mi- 
nus the exciton binding energy {—E^). On the other 
hand the difference between -~E§ and the eigenvalues of 
trion Hamiltonians 17I8|I is equal to the trion binding en- 
ergies {Eg , Eg^) with respect to dissociation into an 
exciton and a free electron (for X^) or a hole (for 
Trion binding energies are equal to the red-shift of the 
trion recombination lines with respect to the exciton line 
in the PL spectrum. 




-3 -2 -1 




12 16 



FIG. 1: (a) [color online] Contour plot of the interaction po- 
tential V = V'mzhi - Zh2) - ^'{L-.ZKi) - V'mzhi) as 
function of the interparticle distances with lateral confinement 
length L = = = 1. Distances and energies are given in 
donor units. The dashed-dotted line corresponds to V = 0. 
(b) The interaction potential plotted for L = 1 along the lines 
Zh2 = —Zhi, Zh2 = —Zhi — 1-5, and Zh2 = —Zhi — 3 marked 
in (a) with (thick) solid, dotted and dashed lines respectively, 
as function of the interelectron {X~) or interhole (X+) dis- 
tance Zi2 = Zhi — Zh2- Thin solid line shows the — 3/|2:i2| 
asymptotic. 

We solve the relative Hamiltonian eigenequations us- 
ing the imaginary time technique2& on a two-dimensional 
grid with a finite-difference approach. We use 201 points 
in both Zhi and Zh2 directions. The size of the computa- 
tional box in both directions is chosen " self-consistently" 
to be 12 times larger than the average distance be- 



tween the particles of the same charge defined as < 

[zhl - Zh2) >^/^- 

In the present calculations we assumed harmonic oscil- 
lator lateral confinement which allows us to simplify the 
problem considerably because of the availability of an- 
alytical formula^'^'^ for the effective one-dimensional in- 
teraction. The confinement lengths le and Ih parametrize 
the strength of the particle localization. Since the single- 
particle energies cancel in the calculation of the trion 
binding energies, the applicability of the present results 
is wider. In fact the present results can be used for any 
form of the lateral confinement (which does not even 
have to be cylindrically symmetric) as long as it produces 
the same effective interaction potential. For instance the 
electron-electron interaction potential for 1^ = 2.95 and 
6 nm is very well (i.e. with a precision better than 2%) 
reproduced for a GaAs quantum wire (me=0.067) with a 
circular square well confinement of depth 320 meV and 
diameters 9.6 nm and 22.8 nm respectively. For elliptical 
harmonic oscillator confinement with different oscillator 
lengths in x and y directions {Ix and ly, respectively) 
we cannot give a closed analytical formula for the effec- 
tive interaction potential. Nevertheless, we have found 
via a numerical integration that the interaction potential 
between two electrons in an elliptical wire can be surpris- 
ingly well reproduced by formula (6) for a circular wire 
with an effective I = {Ix + ly)l^- The numerically calcu- 
lated deviation between the two potentials is not larger 
than 2% for any interelectron distance. The essential as- 
sumption of the present model therefore does not rely on 
the form of the lateral confinement but on its strength, 
which has to be large enough to prevent the Coulomb 
interactions from deforming the lateral wave functions. 
The applied assumption of the frozen lateral degrees of 
freedom for the electron and the hole is applicable for the 
exciton binding energy when Zg < O'd and < a^- This 
condition guarantees that the length of the lateral con- 
finement of the carriers is smaller than the bulk exciton 
radius, and that the sum of the lateral confinement ener- 
gies for the electron and the hole are at least two times 
larger than the exciton binding energy in bulk. For trions 
the applied approximation is better justified and the con- 
ditions are less stringent because the trions have a larger 
size and have smaller binding energies than the exciton. 



III. RESULTS 



We start the presentation of our results by discussing 
the properties of trions in quantum wires with equal lat- 
eral confinement for the electron and the hole (subsection 
III. A), and then in subsection III.B we show the effect of 
different confinements for the stability of X+ and A"^ tri- 
ons. Subsection III.C describes the effect of the magnetic 
field oriented parallel to the axis of the wire. 
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FIG. 2: Wave functions for negative (a,d,g,j) and positive (b,c,e,f,h,i,k,l) trions for le = Ih = L = 0.2 in Zhi and Zh2 coordinates 
(horizontal and vertical axis, respectively) for different values of the mass ratio a. Plots (c,f,i,l) show the wave functions of the 
excited X"*" state antisymmetric with respect to the interchange of the holes. The dashed line in (f,i,l) shows the node of the 
wave functions. Plot (c) corresponds to an unbound state, for other plots the computational box is larger than the fragment 
displayed and the states are bound. 



A. Identical electron and hole lateral confinement 



For equal electron and hole lateral confinement {le = 
Ifi = L) the electron-electron, the hole-hole, and the 
electron-hole interactions have the same form. The total 
interaction potential, identical for both types of trions, 

is plotted in Fig. 1 (a) as fimction of z^i and Zh2 for 
L = 1. The regions of positive (negative) potentials are 



plotted with red (blue) colors. Zero of the interaction 
potential is marked with a dash-dotted line. The inter- 
action potential is minimal along the lines Zhi = and 
Zh2 = at which one of the two electrons and the hole are 
in the same position (for X^) or the position of one of 
the holes coincides with the electron position (for X~^). 
The potential is maximal along the diagonal Zhi = Zh2 at 
which the two particles of the same charge arc localized 
in the same point along the wire length. Fig. 1(b) shows 
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FIG. 3: Electron-hole (solid lines), electron-electron, hole- 
hole (dashed lines) pair correlation functions plots for X~ 
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FIG. 4: Binding energies of the negative (dashed lines) and 
positive trion (solid lines) states for L = 0.2 as function of the 
mass ratio a. Higher solid curve corresponds to the X'^ state 
antisymmetric with respect to the interchange of electrons and 
holes, ie. it is the first excited state. Thin vertical lines show 
the values of a = 1, 1.98, 6.72 and 15.2. Dotted curve, referred 
to the right axis, shows the exciton ground-state eigenvalue. 
Energies and lengths are in donor units. 



the cross-sections of the interaction potential along the 
three straight lines in Fig. 1(a) as function of the inter- 
electron {X~) or interhole (X+) distance zi2. On the 
antidiagonal the interaction has the form of a triangular 
potential well [cf. solid line in Fig. 1(b)]. Moving along 
the antidiagonal is equivalent to interchange the position 
of the two particles of the same charge with fixed posi- 
tion of the third particle of the opposite charge. Along 
the paths plotted with dotted and dashed hues in Fig. 
1(a), which are shifted below the antidiagonal, the po- 
tential has the form of a double potential well [see Fig. 
1(b)] with a barrier near the diagonal resulting from the 
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FIG. 5: [color online] Binding energies of the trions as func- 
tions of the length of the lateral confinement. Lines for 
a — 1, 1.98, 6.72 and 15.2 plotted with black, blue, red, and 
green colors, respectively. For a = 1 binding energies of X~ 
and X^ are equal. The dashed lines are the energies of the 
excited X"*" states antisymmetric with respect to the inter- 
change of the holes {X~ does not possesses a bound excited 
state for a > 1). Lines for X~ at cr = 1.98 and 6.72 have been 
omitted for clarity - they are situated between the a — 1 line 
and a — 15.2 line for X~(see Fig. 01 . Inset: energy eigen- 
values for neutral exciton and charged trions for a — 15.2. 
Energies and lengths are given in donor units. 



repulsion of the equally charged particles. For large zi2 
the potential approaches — 3/zi2 asymptotically which is 
shown by the this solid line in Fig. 1(b). 

Contour plots of the wave function of the negative and 
positive trions calculated for different effective mass ra- 
tios are plotted in Fig. 2 for L = 0.2. For Zg = Ih the 
negative and positive trion relative Hamiltonians (7) and 
(8) differ only by the factor standing in front of the mixed 
derivative (l/cr for and 1 for X+). The ground-state 
wave function for cr = 1 is the same for both trions [cf. 
Figs. 2(a,b)]. The first excited state of X~^, which is an- 
tisymmetric with respect to the interchange of the holes 
is unbound. Its wave function calculated for the size of 
the computational box 100 x 100 (in donor Bohr radius 
units) is plotted in Fig. 2(c). For the unbound state no 
computational box is large enough (i.e. the wave func- 
tion vanishes only at the ends of the computational box). 
The wave function is nonzero only near both axis. One 
of the holes stays at the position of the electron and the 
other strives to be as far as possible from the other two 
particles. 

Results for a — 1.98 plotted in Figs. 2(d-f) cor- 
respond to CdTe material parameters {mh = 0.19mo, 
TOg — 0.096mo) with the donor units ad = 5.4 nm and 
2Rd — 27.6 meV. The probability density maximum for 
(7 = 1 is split into two extrema at the antidiagonal of 
the plots [cf. Figs. 2(a-b)]. For a > 1 these two extrema 
merge into a single one for X~ [see Fig. 2(d)] and for X~^ 
they become more distinctly separated [see Fig. 2(e)] and 
the excited state for X+ becomes bound [cf. Fig. 2(f)]. 
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FIG. 6: [color online] Contour plot of the interaction potential for (a) negative trion V = V^^ {le; Zhi — Zh2) — V^leh', Zhi) ~ 
V°^{leh', Zh2) and (b) positive trion V = V^^ {lh \ Zhi — Zhi) — V^^ {leh \ Zhi) — V^leh', Zh2) as function of the interparticle distances 
for the lateral confinement lengths U = 1 and Ih ~ 0.5. Distances and energies are given in donor units. 
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FIG. 7: [color online] Shifts of the trion recombination PL 
lines with respect to the PL exciton line as function of the 
hole confinement length ilh) for GaAs. Different values of the 
electron confinement length are plotted with different colors. 



Figs. 2(g-i) and 2(j-l) for a = 6.7 and cr = 15.2 
correspond to GaAs (m/i — 0.45too, THe — 0.067mo, 
2i?d = 11.9 meV, — 9.8 nm) and InAs [m-t^ — 0.41mo, 
me = 0.027mo, 2i?d = 3.2 meV and ad = 29.7 nm) ma- 
terial parameters, respectively. Increasing a has an op- 
posite effect on the X~ and wave functions. For the 
negative (positive) trion the local minimum along the di- 
agonal Zhi = Zh2 is less (more) pronounced. The wave 
function evolution with a is related to the tunnelling of 
the particles of the same charge via the potential barrier 
presented in Fig. 1(b). The electrons in with light ef- 
fective masses tunnel easily through the diagonal barrier 
due to the interelectron repulsion. On the other hand the 
diagonal barrier is effectively much larger for the heavy- 
mass holes which prevents its penetratation at large a 
which leads to the appearance of the characteristic max- 



ima elongated along the diagonal in Figs. 2(e,h,k). 

The correlation between the particles in the complexes 
is more clearly visible in the pair correlation functions 
plotted in Fig. 3. The electron-hole correlation function 
is calculated as feh{z) = J dzhidzh2\x{zhi, Zh2)\'^S{z - 
Zhi) and the hole-hole (for X~^) , and electron-electron 

(for X-) as fsanie{z) = ^ dzhldzh2\x{zhl, Zh2)\^ 5[z - 

{zhi — Zh2)]- The Coulomb hole in the hole-hole correla- 
tion in X'^ is much larger than for electrons in X~ but at 
the expense of slightly weaker electron-hole localization. 

The binding energies of the exciton and the trions for 
L = 0.2 are plotted as functions of a in Fig. 4. All the 
binding energies are increasing functions of a. In bulk 
the first excited state of the positive trion is antisym- 
metric with respect to the hole interchange,'^^ possess 
the P symmetry and is bound for a > 4.2. The criti- 
cal value of the mass ratio is much smaller for quasi ID 
confinementi^ Here, for L = 0.2 the excited X+ state is 
bound for cr > 1.2 (see Fig. 3). For quasi ID confinement 
the lowest excited state has the S symmetry with respect 
to the axis of the wire but is of odd spatial parity, i.e., 
it is antisymmetric with respect to simultaneous change 
of sign of all the z coordinates (see Fig. 2). The ground- 
state of Ar+ becomes degenerate with respect to the sym- 
metry of the wave function, i.e., the hole interchange, for 
large a for which tunnelling through the diagonal poten- 
tial barrier (cf. Fig. 1) disappears. At large a also the 
probability density of the excited AT"*" level becomes iden- 
tical to the ground-state probability density (cf. Fig. 2). 
X~ does not possesses a bound excited state for cr > 1. 

The inset to Fig. 5 shows the ground-state energy for 
cr — 15.2 as function of the lateral confinement length. 
In the L — limit the average interparticle distances de- 
crease to zero and the energies diverge to minus infinity. 
This is a consequence of the Coulomb interaction singu- 
larity in one dimension. ™ The main part of Fig. 5 shows 
the shifts of the trion PL lines with respect to the exci- 
ton line (calculated as the difference of the eigenvalues 
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presented in the inset) for different values of cr. It turns 
out that the binding energies have a power law depen- 
dence on L, i.e. i^"^, for the and ground state 
presented in this figure q changes from 0.83 (cr = 1) to 
0.91 {X+ for a = 15.2). 



B. Effect of different electron and hole lateral 
confinement 

Let us now consider the interaction potential for 
stronger hole confinement. Figs. 6(a) and (b) show the 
interaction potentials for = 1 as in Fig. 1 but for 
smaller Ih = 0.5. For both the negative [cf. Fig. 6(a)] 
and the positive trion [cf. Fig. 6(b)] the potential min- 
ima at Zhi = and Zh2 = become deeper with respect 
to the le = Ih case presented in Fig. 1. For X~ the 
electron-electron interaction (the diagonal potential bar- 
rier) is not affected by the change of lu [compare Fig. 1 
and Fig. 6(a)]. On the other hand the hole-hole repulsive 
interaction for X^ is strongly increased. 




/g [nm] 

FIG. 8: [color online] Difference of the positive and negative 
trion binding energies (in meV) as function of the electron 
and hole confinement lengths for GaAs material parameters. 
Blue (red) regions correspond to more stable negative (posi- 
tive) trion. Above the dashed-dotted line the negative trion is 
unbound. The green line corresponds to Eb{X~) = 4.2 meV 
and the yellow line to Eb{X'^) — 2.9 meV. 

The effect of the hole localization on the trion binding 
energies is plotted in Fig. 7 for GaAs material parame- 
ters and fixed values of the electron lateral confinement. 
Consistently with the results of section III. A for l,, — Ih 
the positive trion is more stable than the negative trion. 
A decrease of In below the value of le results in the in- 
terchange of the X~ and X+ energy hnesiSS This is due 
to the enhanced hole-hole interaction shown in Fig. 6 
(b). The negative trion binding energy is a monotonous 
function of the hole confinement length, the larger 1^ the 
smaller is the electron- hole interaction stabilizing X^ . 
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FIG. 9: Interparticle distances for the negative (solid lines) 
and positive trions (dashed lines) and for the exciton (dotted 
line) for GaAs material parameters and le ~ 2.95 nm. (e-e), 
(e-h), and (h-h) stand for the electron-electron, the electron- 
hole and the hole-hole distance. 

The situation is more complex for X'^, since with in- 
creasing Ih also the destabilizing hole-hole interaction 
decreases. As a consequence the positive trion binding 
energy possesses a maximum as function of Ih- 

The difference of the positive and negative trion bind- 
ing energies is plotted in Fig. 8. Both the trions are 
equally stable for Ih = 0.92Ze — 0.38 nm. For Ih larger 
(smaller) than 0.921^ — 0.38 nm X~^ is more (less) sta- 
ble than X^ . The maximum of the AT^ binding energy 
presented in Fig. 7 follows a path that is nearly lin- 
ear for le > 2 nm and is approximately parametrized 
by Ih ~ 1.62le — 1.98 nm. For the points at the left of 
the dash-dotted line the electron-hole attractive interac- 
tion stabilizing the trion is so weak with respect to the 
electron-electron repulsive interaction destabilizing the 
complex that the negative trion stops to be bound (see 
also the hue for X~ at le = 2.95 nm in Fig. 7). The 
absence of the negative trion binding requires a substan- 
tially weaker hole localization than the localization of 
the electron which is rather impossible to obtain in the 
presently produced quantum wires and would require the 
valence band offset between the wire and the matrix to 
be much smaller than the conduction band offset. More- 
over, the present modelling based on the assumption that 
the lateral wave functions are not affected by the inter- 
action is likely to fail since the hole wave function is very 
likely to become more localized due to the attraction by 
strongly confined electrons. 

The fit of the calculated X^ and X^ binding energies 
to the experimental data is obtained at the crossing of 
the green and yellow lines, i.e., for le — 2.95 nm and Ih — 
1.3 nm. The obtained fit corresponds to realistic values 
which give a general idea on the particle localization in 
the wire (the measurementsii were performed on a V- 
groove GaAs/AlGaAs quantum wire with a thickness of 
the GaAs crescent of 3 nm at the center). Obviously, a 
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more realistic model is required to extract details of the 
confinement from the experimental data. 

The dependence of the size of the trion, i.e., the in- 
terparticle distance as function of the hole confinement 
length, is shown in Fig. 9 for 1^ = 2.95 nm. The electron- 
hole distance for the trions and the exciton have been cal- 
culated as ■y/<z^j>, and \/<z^h^, respectively. The 
hole-hole distance for X+ and th e electron-electron dis- 
tance for X~ are determined as ^ < {zhi — >■ The 
size of the exciton increases as 1^ increases which is due to 
the reduced value of the electron-hole interaction. Much 
stronger dependence on Ih is observed for , which be- 
comes unbound for Ih > 7 nm [cf. Fig. 7]. The de- 
pendence of the X'^ size on the hole confinement is non- 
monotonous. The positively charged complex has the 
smallest size near Ih — S nm when it is the most strongly 
bound [cf. Fig. 7]. For Ih — le ^ 2.95 nm the order 
of the interparticle distances in the two complexes is the 
same as in two-dimensional quantum wells (compare Fig. 
4 of Ref. fi2i] ) . In spite of the fact that the probability 
of finding both holes in X"*" in the same position is much 
smaller than for electrons in X~ (cf. Fig. 3) the longer 
tail of the electron-electron correlation function results 
in a larger electron-electron distance than hole-hole dis- 
tance. 




B[T] 



length {Ic = 36.28/^/B nm VT). Since 1(3) for high 
magnetic fields decreases to Ic, the magnetic field tends 
to equalize the electron and hole localization. The bind- 
ing energy of the trions can be obtained following paths 
of {lh{B),le{B)) on Fig. (8). The modification of the 
binding energies of the trions by the magnetic field is 
presented in Fig. 10 for different oscillator lengths. In a 
magnetic field of 40 T, 1{B ~ Q) ~ 4.9 nm is decreased 
to 1{B = 40T) = 4.4 nm and 1{B = 0) = 6.86 nm to 
1{B = 40T) = 5.19 nm. For 1^ = Ih = 6.86 nm (see the 
red curves in Fig. 10) the magnetic field decreases the 
length of confinement exactly as in the case presented in 
Fig. 5. In the more realistic case of stronger hole confine- 
ment, i.e., for le = 6.86 nm and Ih = 4.9 nm (see the black 
lines in Fig. 10) the magnetic field increases the binding 
energy of the positive trion more strongly. This is be- 
cause the magnetic field more strongly affects the larger 
le value than the smaller hole localization length Ih which 
increases the electron-hole interaction more strongly than 
the repulsive hole- hole potential. For X~ the effect of the 
increased electron-hole interaction is nearly cancelled by 
the increase of the electron-electron potential value. Near 
35 T a crossing of the black lines is observed, which cor- 
responds to passing from the "blue" to the "red" region 
on the phase diagram of Fig. 8. This crossing is qualita- 
tively opposite to the one obtained for two-dimensional 
quantum wells, ^"^ in which a small magnetic field (around 
1 T) increased the X^ stability over the X+ binding 
energy. In quantum wells the magnetic-field related in- 
crease of the single-particle energy is smaller for electrons 
and holes bound in the trion complex than for the low- 
est Landau level in the final state of the free electron 
and hole after the trion dissociation. The crossing of the 
binding energies observed in quantum wirea^ is there- 
fore at least partially due to the stronger dependence of 
the electron lowest Landau level - note an almost linear 
magnetic field dependence of the trion binding energies 
in Fig. 2 of Ref. [^^]. In the present calculations the free 
electron and the free hole are strongly localized in the 
plane perpendicular to the field and the single-particle 
magnetic field effects cancel due to the assumption of 
the frozen-lateral degrees of freedom, so that the cross- 
ing is entirely due to the modified effective inter-particle 
interactions. 



FIG. 10: [color online] Magnetic field dependence of the trion 
binding energies in GaAs for different values of the electron 
and hole oscillator lengths. 



C. Magnetic field parallel to the wire 

In the present approach it is straightforward to in- 
clude a magnetic field oriented parallel to the axis of the 
wire. It simply scales down the electron and hole oscil- 
lator lengths according to the formula 1{B) = {l/l'^{0) + 
l/lf)-^/^, where Ic = ^/2h/eB is the magnetic field 



In the case of stronger electron confinement {1^ = 4.9 
nm, Ih = 6.86 nm - the blue lines in Fig. 10) the situation 
is just opposite, the X~ is less strongly bound and the 
magnetic field acts more strongly on the negative trion 
decreasing the binding energy difference according to the 
mechanism described above. However, for the electron 
confinement equal or stronger than the hole confinement, 
the magnetic field does not lead to crossing of the trion 
energy lines and X~^ is more stable for any magnetic field. 
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IV. SUMMARY AND CONCLUSIONS 

We studied the properties of the negative and posi- 
tive trions in quantum wires with strong lateral confine- 
ment using the approximation of the lowest subband oc- 
cupancy which allows for a numerically exact solution of 
the multi-particle Schrodinger equation. We investigated 
the relative stability of the positive and negative trions 
with respect to the dissociation into an exciton and a free 
carrier for different electron and hole confinement. We 
found that the order of the negative and positive trion PL 
lines is interchanged when the lateral confinement of the 
hole is stronger than the one for the electron. In a GaAs 
quantum wire with If, = 5 nm we predict that when 
is 20% smaller the positive and negative trion recombi- 
nation lines interchange. The change in the order is due 
to modification of the effective interactions in the trion 
complexes. The present results provide an explanation 



for the recently experimentally observed larger stability 
of the negative trion in quantum wires We predict that 
for larger X~ stability due to stronger hole confinement 
the magnetic field oriented parallel to the axis will tend 
to interchange the order of the X'^ and X~ energy lines. 
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